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Abstract 

The first discrete Painleve equation (dPI), which appears in a model of quan- 
tum gravity, is an integrable nonlinear nonautonomous difference equation 
which yields the well known first Painleve equation (PI) in a continuum limit. 
The asymptotic study of its solutions as the discrete time-step n — > oo is im- 
portant both for physical application and for checking the accuracy of its 
role as a numerical discretization of PI. Here we show that the asymptotic 
analysis carried out by Boutroux (1913) for PI as its independent variable 
approaches infinity can also be achieved for dPI as its discrete independent 
variable approaches the same limit. 
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1 Introduction 



Our aim is to study the equation 

an + (3 . . 

Vn+l + Vn + Vn-i = h 7 (1) 

Vn 

in the limit as n — > oo. Eqn ([!]) is known as the discrete first Painleve 
equation (dPI) because the scaling limit y n = 1 + h(— 2w(x) + c/2), x = 
Vhn - 6/2 - c 2 /24, with a = h 3 / 2 a, f3 = -7 + 3 + h 2 b, 7 = 6 + he, yields the 
classical first Painleve equation (PI) 

w" = 6w 2 + x, 

as h — >• 0. 

PI is the simplest of six well known nonlinear second-order ordinary differ- 
ential equations (ODEs) in the complex plane called the Painleve equations. 
Their characteristic property that all movable singularities of all solutions are 
poles is called the Painleve property. Painleve |J, Gambier p|, and Fuchs 
[0] identified them (under some mild conditions) as the only such equations 
with the Painleve property whose general solutions are new transcendental 
functions. 

More recently, it was discovered by Ablowtiz, Ramani, and Segur j| that 
such equations are closely related to integrable partial differential equations 
[||, Q . There is an arsenal of rich techniques that exhibit the special character 
of the Painleve equations as integrable equations. In particular, Painleve [0 
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showed that they can be written in terms ratios of entire functions, Fuchs 
[||] and later Jimbo et. al || (see also 0) showed they are isomonodromy 
conditions for associated linear ODEs, and Joshi and Kruskal [|K| showed 
that their connection problems can be solved directly. Consequently, their 
solutions play a distinguished role as nonlinear special functions. A natural 
development is to ask: What are their corresponding discrete versions? 
The first insight into this question came from a discrete version of PI de- 



rived in a model of quantum gravity ITT], [L2| . Subsequently, a tour-de-force of 



work was accomplished by Grammaticos, Ramani and collaborators (see [13 



[PH and references in fll5|| ) to answer this question. They derived integrable 
discrete versions of the first to the fifth Painleve equations and identified a 
property, called the singularity confinement property, that appears to be a 
discrete version of the Painleve property. Isomonodromy problems have been 
found for these discrete Painleve equations. They appear, therefore, to be 
integrable and deserve further study as rare examples of nontrivial integrable 
systems. An asymptotic description of their possible orbits as n — > oo is a 
natural part of this study. 

The limit we study, i.e. n — > oo, is related to the continuum limit. In 
such a limit, n is typically multiplied by a small parameter — > 0, where 
\x = 1/2 for dPI to PI, while nh^ remains 0(1). Such a limit can only be 
achieved if n is itself approaching infinity. In order to understand how the 
orbits of dPI relate to the solutions of PI in a continuum limit, we need to 
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study dPI in the limit n — > oo. This is crucial also for deciding whether or 
not dPI is an accurate numerical discretization of PI. 

It should be noted, however, that the limit n —>■ oo is not equivalent to 
a continuum limit. Continuum limits of dPI involve scaling limits of y n and 
the parameters a, (3, 7, in addition to requiring n — > 00. Moreover, there 
exist limits of dPI other than the one described above to PI. For example, an 
alternative continuum limit to PII is also known. Nevertheless, for simplicity, 
we restrict the comparison of our discrete results to those known for PI. 

We show that the asymptotic behaviours admitted by dPI as n — > 00 
are qualitatively extremely similar to those found for PI by Boutroux [|i~6[| . 
Boutroux showed that in general the asymptotic behaviours of the (first four) 
Painleve transcendents near infinity in the plane of the independent variable 
are given by elliptic functions. He studied the behaviour of the transcendents 
along a line of approximate periods of such elliptic functions and he showed 
that within certain sectors near infinity there also exist special behaviours 
which are asymptotically free of poles. His deductions were facilitated by 
a transformation of variables that made the asymptotic behaviours more 
explicit. For PI, Boutroux's transformation is 

w(x) = \fxu{z) 

z = -x 5 ^ (2) 
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under which PI becomes 



where the primes now denote differentiation with respect to z. (See JLO 
for a deduction of Boutroux transformations based on maximal dominant 
balances.) It is clear that on a path where u, u', u" are bounded, and u" is 
of order unity, the leading-order behaviour of Eqn(|3|) as \z\ —* oo is given by 

u" ss 6u 2 + I, 

which is solved by Weierstrass elliptic functions. On the other hand, if u" <C 
1, as \z\ — > oo, the leading-order behaviour is given by 

2 1 
u ~ , 

6' 

which gives rise to pole-free behaviours. 

There exists a Boutroux transformation of dPI that also facilitates its 
asymptotic study as n —>■ oo. An analysis of possible dominant balances of 
dPI as n — > oo shows that there is only one maximal balance 

an 

Vn+l + Vn + Vn-1 ~ , 

Vn 

which is achieved when y n grows like 0(\/n). Hence we are led to 

Vn = Vnu n 



5 



which transforms dPI to 

a 7 (3 

u n +i + u n + u n -i = h — == H 

u n Vn nu n 

fy/n+1- y/n\ j \Jn - 1 - y/n \ 

~ { — ^ — S Un+1 ~{ — — S Un - UA) 

In Section 2, we show that this equation is solved to leading-order by elliptic 
functions generally and by pole-free behaviours less generally. 

The elliptic-function-type behaviours are described by an energy-like pa- 
rameter E and a phase-like parameter $ which are constant to leading-order 
but vary slowly as n — > oo. The evolution of E over many periods of the 
leading-order elliptic function can be studied via an averaging method. We 
carry out this study in Section 3. 

Although qualitatively we obtain very similar results to those found by 
Boutroux, we point out here that there are differences. In the special case 
7 = 0, the main result of Section 3 can be regarded as the discrete analogue 
of the corresponding result found by Boutroux for PI. However, in the case 
7 7^ 0, the slow evolution of E has an additional component. This additional 
complexity is due to the fact that our asymptotic limit n — ■> oo is not exactly 
a continuum limit. Moreover, the limit n — > oo includes the case where 
x = nh^ is finite. Our analysis therefore includes regions where Boutroux's 
asymptotic analysis may not be valid (he assumed x — > oo). 

What we demonstrate in this paper is that a comprehensive (formal) lo- 
cal asymptotic description of dPI can be carried out as n — > oo. Although 
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our results are obtained as n — > oo, a similar leading-order analysis can be 
carried out in neighbourhoods of any ordinary point n = uq. This also yields 
elliptic-function-type behaviours and pole-free behaviours. Hence our analy- 
sis provides futher evidence of the special character of dPI as an integrable 
equation. Ramani, Grammaticos et al have pointed out that dPI has the sin- 
gularity confinement property by a perturbation analysis carried out around 
its singularity at y = 0. Our asymptotic results verify this pole-like singular- 
ity structure of the orbits to leading order by showing that locally in patches 
near infinity, the orbits are given by meromorphic functions. 

2 Local Asymptotic Results 

Here we analyse Equation (|J) to find the locally valid asymptotic behaviours 
of dPI in the limit n — > oo. Our results are qualitatively the same as those 
found by Boutroux for PI. 

Suppose that u n , u n+ i, u n _i are all 0(1). Then we have a maximal 
dominant balance of Eqn (|4]) given by 

U n +1 + U n + U n -i « . (5) 

Henceforth we adopt the notation 

u n+1 =: % u n = u, u n _i =: u. 
Multiplying Eqn by uu and uu respectively and subtracting the two 
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results, we get 



uu 2 + u 2 u — u 2 u — uv 2 



a(u — u 



Both sides of this equation are exact differences and so we can integrate 
(actually sum) to get 



Note that E is constant to leading-order as n — > oo, but varies slowly with n. 
That is, E — E is not zero but small (in fact, of order 0(1/ y/n) for n ^> 1, see 
Eqn (14)). We will refer to E as an energy-like integral (or parameter) for 
dPI. Consider the left side of Eqn(|6]) as a polynomial in x = u, y = u. The 
value of E (to leading-order) defines level curves of this polynomial which 
are parametrized by elliptic functions (see [0). 

An alternative way to see that the solution involves elliptic functions is 
to sum the equation once more. Solve Eqn ([]) for u and write the result as 



uu 2 + u 2 u — a(u + u) =: E. 



(6) 




where 



P(u,E) := (u 2 + a) 2 + AEu. 



Then the result can be written as 



u(u - u) 



1 



which upon integration gives 



E 



rt 



u k (u k+ i - u k ) 



n — $ 



k=l 
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where $ is a constant. The left side is a discrete analogue (Riemann sum) 
of an elliptic integral of the first kind. 

The leading-order elliptic function has two periods. The periods are (in 
general) complex numbers given by the fact that the curve defined by Eqn 
(H) has genus two. In other words, P has four branch points (for generic 
values of E) and there exist two linearly independent closed contours that 
enclose a pair of these. We will denote the two periods (given by the two 
contours Cj) as Uj, j = 1, 2. 

In fact, the elliptic function is a discrete sampling of the continuum (Ja- 
cobian) elliptic functions. So the definition of Uj requires slightly more ex- 
planation. We regard them as being defined by the elliptic integral 



where the path of integration for each given E is a closed contour (as de- 
scribed above) formed by interpolation through the orbits of u parametrized 
by the phase 

The elliptic-function-type behaviour u of dPI, possesses these periods only 
to leading-order. However, we nevertheless define u}j(E) in the same way and 
regard them as (implicit) functions of n through their dependence on E. We 
will also call them (loosely) periods. 

There exist special values of E for which the elliptic functions degenerate 
to singly periodic (trigonometric) functions. To find these values, rewrite P 
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as 

\2/2 



P(u, E) = (u- p) 2 (u 2 + au + r). 
Multiplying out the product on the right and equating coefficients gives 

a = 2p , r - 2pa + p 2 = 2a (8) 
p 2 a - 2pr = 4E , p 2 r = a 2 (9) 

Solving these shows that 

p 2 = a/3, —a 

and 



E = —a y/a/3, 

respectively. For these values of E, the elliptic integrals described above 
degenerate to yield only a single period for u. 

Now suppose that u, u, u are still 0(1) but u ~ u, u ~ u. Then the 
first integral (^|), although still defined, describes an algebraic function (to 
leading order). Since u is constant to this order, where the constant is given 
by 

a 

3u w -, 10 

u 



we write u = y ( a /3) + f , where t)<l. Then Eqn (^) becomes 

u + 4u + ?; = -t=(l + v/J(a/3)) 

+0(^l/n,vv } v 2 ,vv j v/n } v/rij (11) 
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We get the asymptotic solution 

^/W3) + ^= + 0(l/ n ). (12) 

This solution obviously has an infinite algebraic asymptotic series expansion 
and, therefore, has no poles for n ^> 1. For that reason, we call such solutions 
pole-free behaviours. 

A straightforward calculation shows that a perturbation of a solution V 
with such behaviour, i.e. u — V + v, gives 

where 

A± = -2 ± V3, 

and c± are arbitrary constants. Obviously, |A_| > 1 and |A+| < 1. So to get 
a consistent perturbation, we must have c_ = 0. In other words, we get an 
exponentially small perturbation, multiplied by c+, which is hidden beyond 
all orders of the algebraic expansion (|l"2"D as n — > oo. Since the free parameter 
c + cannot be identified uniquely, the pole-free behaviour ( |T2"D cannot describe 
a solution uniquely. 

We chose the solution given by A + above under the assumption that n is 
real, positive. However, in general, dPI could be regarded as a mapping posed 
along a complex line, i.e. n could be complex. Therefore, the perturbation v 
could be given by either (A±) n or indeed a combination of both (e.g. if n is 
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pure imaginary). In other words, the behaviours given to the first few orders 
by ( |T2"D suffer from Stokes' phenomenon in the complex n-plane. 

There is one other dominant balance giving rise to pole-free behaviours. 
In obtaining V above, we assumed that 

u ~ c, u + u P3 2c, 

and then found the constant c to be given by J (a/3). An alternative con- 
sistent assumption is 

u ~ c, u + m ~ —2c. 

In this case, we find that 

u « (±l) n v /z a + 0(l/ v / ^). (13) 

Again the solutions with this leading order behaviour suffer from Stokes' 
phenomenon. However, now the Stokes' lines are orthogonal to those of the 
pole-free behaviour given by V. Writing u = W + w, where W is given to 
leading order by Eqn([L3|), we get 

w ss b±(±i) n , 

where b± are arbitrary constants. So along purely real directions in the n- 
plane, the perturbation w is approximately oscillatory, whereas along purely 
imaginary directions, one of b± must be zero and w is hidden beyond all 
orders of the divergent expansion for W. 
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Along anti-Stokes' lines, i.e. directions along which the perturbations v 
or w are approximately oscillatory, the perturbations are not simply linear 
combinations of the exponentials found above. As for PI (see |L8|, |10ll), they 
are more accurately represented by (discrete) Fourier series which yield an 
expansion composed of increasing powers of such exponentials multiplied by 
algebraic prefactors. We omit the details here for simplicity. 

It can be shown that the elliptic-function-type behaviours and the pole- 
free behaviours are connected in the space of solutions. The latter can be 
shown to be degenerate limits of the former, attained when 



Although the number of pole-free behaviours of PI (or PII) differs from 
that for dPI (PI has two, PII has three), these results reflect qualitatively 
exactly the asymptotic results found by Boutroux. 

3 Averaging 

In this section, we analyse the change of the energy-like parameter E (describ- 
ing the elliptic-function-type behaviours) over a period ujj of the leading-order 
elliptic function by using an averaging method. 

E fluctuates with n. This is clear from the definition @ which gives 




E-E 




(14) 
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Suppose we are given initial values at some point |$| 3> 1, which define 
a solution as follows: 



Assume the periods O0j to be those given at <£> by this value of E. 

To carry out an averaging method, we assume that there is a smooth, 
slowly varying function interpolated between the successive points on the 
orbit defined by E and $ (or the initial values above). We wish to study the 
slow variation of E from one period of the leading-order elliptic function to 
the next. Let u be 



where U represents the slowly-varying part of u and s -C U represents its 
fast fluctuations (as |n| — > oo). The initial values are now 



u($) = l, u($ + l)=p(^0). 



Then we have 



J5($) =p 2 + p - a (p+ 1). 



u = U + s, 



1 







s($ + 1) = 0. 



(15) 



Then Eqn (Q) gives 



u + u + u 



a 



(16) 
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-(^l + l/n-l)(U + s)-^/l-l/n-l\(U + s) 

(17) 



Note that U is an elliptic function (by arguments presented in Section 2). 

Consider the first integral given by Eqn(||). Expanding the solution u as 
U + s and keeping only terms of order s or 1/ y/n gives 

s(2UU + U 2 -a) + s(JJ 2 + 2UU - a) = 1 ^=. (18) 



This equation can be integrated once more by using the following observa- 
tions. The defining equation for U gives 

UU-UU = UU + U(U + U -a/U) (19) 

= 2UU + U 2 -a (20) 

-VV+UU = U(U + U - a/U) + UU (21) 

= U 2 + 2UU-a (22) 



Letting 



we have 



Hence Eqn([18D becomes 



F := UU — UU, 



F = UU -UU. 



^ ^ UU 
sF -sF = 7- 



\/n' 
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There is an integrating factor for this equation. After division by FF, it 
becomes 



s 
F 



1-- 



UU 



FFy/n 

Hence we can integrate (actually sum) to get s as 



_ jp v> Uk+iUk 



fc=$ F k+ iF k \/~k 

where the sum is understood to be zero in the case n — 3>. 

Now consider EqnflT^D for E. Note that the first term on the right gives 

UU 



E(<S> + lu) - E($) = 7 



-0{l/$,s/V$), 



to leading order. However, the contribution of this term is actually smaller 
than 0(1/ y/n) because U is periodic and 

1 1 



O (l/$ 3/2 . 



To get the nonzero slow change of E, therefore, we need to consider terms of 
size 0(l/n, s/\/n). 



At these orders, we get 
E($ + lu) -E($) 



7- 



Us + Us 



(3{U + U) 



n 



•4 £ u k {u k+x -u k _ x f 
/v fc=<i>+i 



(23) 
(24) 
(25) 
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By periodicity and integration by parts, the second ratio on the right can be 
shown to be zero to leading order. Hence we have 



E($ + u) - E($) 



7 2 [ TT ^ ^4 U k+1 U k x ^ ( /,+ [( /, 



Uk+lUk 



1 



— 53 u k (u k+1 - u k ^ 



2$ 



$+1 



fc=$ 



1 

2$ 



(26) 



$+1 



To simplify this expression, we use the notation P(u, E) defined in Section 
2. Note that 



F = U(U -U) = JP(U, E) 



Therefore, we have 



UU 

TP 



UU 



P(U,E)P{U,E) 



Hence Eqn(p6|) becomes 



E($ + uj) 



Uis> + iF<s, + C/$F$ + i ) 



U k+1 U k 



$ JP(U k+1 ,E)P(U k ,E) 



-— y 

$+1 



P(U k ,E) 



2$ ^ f4 



In the case 7 = 0, this result yields 
£($+cj)-£($) «- 



1 *+?P(U k ,E) 



E 



2$ frJ I/* 



(27) 



(28) 
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This is the discrete analogue of the result found by Boutroux. To recognize 
this, note that 



and so the sum on the right of Eqn (|28|) is a discrete analogue of the elliptic 
integral 



(Note that the derivative of this integral with respective to E is the previously 
described elliptic integral ujj and that in making continuum analogies we have 
ignored numeric factors.) This is what Boutroux found (see page 319 of fT6H ) 
for the evolution of E (called D in his paper). He showed by using this result 
that along a chain of increasing points $j defined by $o = Eq — 
$j = + u(Ej-i), Ej = E($j), E must be bounded and approach a value 
for which the leading order elliptic function degenerates to a trigonometric 
function. 

The generic case 7 ^ differs markedly from this (see the remarks made 
near the end of Section 1) and will form the subject of future explorations. 
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P(U k ,E) 



(U k+1 -U k ^)y/P(U k ,E), 
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